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Abstract 

In the present work we examined the hypothesis that, a core mass func- 
tion (CMF), such as the one deduced for cores in the Orion molecular cloud 
(OMC), could possibly be the primogenitor of the stellar initial mass function 
(IMF). Using the rate of accretion of a protostar from its natal core as a free 
parameter, we demonstrate its quintessential role in determining the shape of 
the IMF. By varying the rate of accretion, we show that a stellar mass distri- 
bution similar to the universal IMF could possibly be generated starting from 
either a typical CMF such as the one for the OMC, or a uniform distribution 
of prestellar core masses which leads us to suggest, the apparent similarity in 
shapes of the CMF and the IMF is perhaps, only incidental. The apodosis of 
the argument being, complex physical processes leading to stellar birth are 
crucial in determining the final stellar masses, and consequently, the shape of 
stellar mass distribution. This work entails partial Monte-Carlo treatment of 
the problem, and starting with a randomly picked sample of cores, and on the 
basis of classical arguments which include protostellar feedback and cooling 
due to emission from warm dust, a theoretical distribution of stellar masses 
is derived for five realisations of the problem; the magnetic field, though, has 
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been left out of this exercise. 
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1. Introduction 

While the formation of stars is now reasonably well understood, the ques- 
tion related to the distribution of stellar masses at birth characterised by the 
so called stellar initial mass function (IMF), and in particular, about the 
shape of the IMF, is still far from being resolved. The multiplicity of stars is 
also well known, in fact, nearly half the stars in the solar neighbourhood have 
at least one companion, of which, the M-dwarfs form a large proportion (e.g. 
Duquennoy & Mayor 1991; Fischer & Marcy 1992; Thies & Kroupa 2007). 
An empirical power-law distribution for stellar masses was first suggested 
by Salpeter (1955), however, this distribution was defined for relatively high 
mass (greater then 1 M0) stars only. More recent observations have revealed 
stellar objects with sub-solar masses, as low as as ~0.02 Mq, the Brown- 
dwarfs (BDs), which also appear to have companions. In fact, the IMF for 
BDs suggested by Thies & Kroupa (2008) shows a discontinuity at the hy- 
drogen burning limit. 

Numerous authors have observationally derived the IMF for massive star- 
clusters such as the IC348 (Lada & Lada 1995; Luhman et al. 1998; Muench 
et al. 2003), NGC 1333 (Aspin, Sandell & Russell 1994; Lada et. al. 1996), 
and Ophiuchus (e.g. Bontemps et al. 2001), to name a few. The similarity in 
shapes of these respective system IMFs has encouraged the idea of a universal 
IMF. A piece- wise form of such an IMF was proposed by Kroupa (2002), and 



2 



a lognormal IMF, by Chabrier (2003). The question relating to the form of 
the IMF raises a more fundamental issue, that of the evolution of prestellar 
cores, and the factors likely to affect the final distribution of stellar masses. 
Infrared observations of dense prestellar cores have enabled determination 
of core masses in a number of nearby star-forming regions, which like the 
stars themselves, also follow a power-law distribution (e.g. Motte, Andre & 
Neri 1998; Andre et al. 2004; Alves et al. 2007; Nutter & Ward-Thompson 
2007, Enoch et al. 2008; Rathborne et al. 2009). The similarity between the 
distribution of core masses and stellar masses is suggestive of the possibility 
that the stellar IMF may derive its form from the former distribution, the 
so called core mass function (CMF) (e.g. Motte et. al. 1998, Nutter & 
Ward-Thomson 2007). 

It has recently been shown by Swift & Williams (2008) that the IMF is 
fairly robust to any variations in the evolutionary history of the prestellar 
cores. Similarly Goodwin et al. (2008), starting from a power-law CMF, 
derived the stellar IMF by assuming a core-to-star efficiency of about 15 %. 
This problem has also been examined numerically by for instance, Bonnell 
et. al. (2006), Clark et. al. (2008), who have shown that gravitational 
fragmentation of prestellar clumps could produce a stellar mass distribution 
similar to the Salpeter IMF, with a turnover located at approximately the 
Jeans mass of the precollapse clump, and a knee at ~0.6 M©. The numerically 
derived IMF by the above authors was shown to become shallower below its 
knee, consistent with the universal IMF suggested by Kroupa. 

In the present work we use a sample of cores with randomly picked masses 
to derive the stellar-mass distribution for five realisations of the problem. 
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First, starting with a power-law CMF defined by Eqn. (1) below, we derive 
the stellar mass distribution for 4 cases: (a) cores in this case are initially 
quiescent, self-gravitating cores arc allowed to cool according to a suitable 
cooling-law, (b) similar to (a), however, precollapse cores are additionally 
supported by a transonic turbulent field, and (c) as in (a), the cores are 
quiescent, but self-gravitating cores are not allowed to cool. Finally, (d) for 
a sample of cores having a uniform mass distribution, we derive a stellar mass 
distribution for quiescent cores while admitting a cooling scheme. This case 
is repeated for a higher rate of mass accretion. The scheme employed for the 
purpose is outlined in §2 below; the results are presented and discussed in §3 
and §4, and we conclude in §5. 

2. Generating the stellcir Initial mass function 

We wish to test how crucial is the nature of the initial distribution of 
prestellar cores in determining the shape of the distribution of stellar masses, 
in particular the resemblance of the latter with the universal IMF. This 
hypothesis is tested using an analytic scheme based on the following three 
premises : 

1. Assumption of the validity of Larson's scaling relation defined by Eqn. 
(2) below, down to the spatial extent of a prestellar core. Keeping in 
mind the smallest length-scale over which the relation is valid, (typi- 
cally 0.1 pc), our sample of prestellar cores is such that the size of the 
smallest core in it is larger than this minimum length scale. 

2. We do not include contribution due to compressional heating of the 
gas in a self-gravitating core, nor do we account for viscous heating of 
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the protostellar disk. Any additional heating will only raise the mass 
of putative fragments and probably produce fewer fragments, causing 
the final distribution to move towards the high-mass end. However, 
it is not our aim in a demonstrative calculation such as this one, to 
resolve the population of sub-stellar dwarfs, a largely recondite area in 
the theory of star-formation. 
3. Investigation of the temporal stability of the attendant circumstellar 
disk is beyond the scope of the proposed deduction, so we simply work 
out the maximum number of sub-stellar fragments likely to condense 
out of the disk over the time required for the protostellar envelope to be 
blown off via feedback. Having accounted for protostellar feedback from 
the first generation of protostars, the calculations were terminated on 
formation of the next generation of stars, where possible, out of residual 
gas in a core. It is reasonable to so do since most of the gas in a core 
will have been consumed by the two generations of stars; some gas from 
the core may indeed be lost via feedback. 

We note that under the proposed scheme, discussed in §2.1, we are only 
estimating the maximum number of fragments likely to condense out of a 
self-gravitating core. In the first three realisations of the problem, we use a 

sample of 14,000 randomly picked cores having individual mass, Mcore/MQ, 
in the range 0.5 to 10. The core masses are distributed according to a power- 
law. 




(1) 



5 



0.5 



2 



4 



8 



Mcore/[Ms„nl 

Figure 1: The histogram plotted on a logarithmic scale shows the initial distribution of 
prestellar cores, according to the power-law defined by Eqn. (1) below. 

with slope, 



suggested by Nutter & Ward- Thompson (2007), for prestellar cores in the 
Orion Molecular cloud; dN is the number of cores in the range {Mcore: Mcore+ 
dMcore)- The initial distribution of core masses in our test sample is sliown 
in Fig. 1 above. For the purpose we used a random number seed, 7^(0, 1), 
that produced a uniform distribution of random numbers between and 1. 
The sample of cores within the desired range of masses is then generated 
by appropriately scaling the uniform distribution, TZ. The radius of a core, 
Rcore, is calculated using the empirical length scaling relation due to Larson 
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Under the approximation of a uniform density sphere for a core, the ensemble 
of prestellar cores in our sample has initial density typically on the order of 
a few times 10^ cm~^, modulo small variations by a factor of 1-2. The mass 
of putative fragments, Mf^ag, the number fragments, and the efficiency of 
fragmentation, i.e. the star- forming efficiency (SFE), were calculated for 
twelve choices of core temperature, Tcore, in the range 7 K to 30 K. Then, 
in a second realisation of the problem, we raised the effective sound-speed 
within a precollapse core by superposing a transonic turbulent velocity field 
having a mean dispersion, a, and given by the scaling relation 



While acknowledging the revision to this scaling relation recently proposed 
by Ballesteros-Paredes et al. (2010), we are of the opinion that it would 
largely serve to increase the magnitude of a. The canonical prescription due 
to Larson, in our opinion therefore, is sufficient for the proposed demonstra- 
tive exploration. 

Fragmentation of a prestellar core The problem of a self-gravitating 
core has been studied analytically with great rigour for close to four decades, 
beginning the seminal work by Larson (1969), Penston (1969), followed by 
Hunter (1977), Shu (1979), Whitworth et al. (1998), and Krumholz et 
a/. (2005). It has also been studied numerically in the recent past by Bonnell 
et a/. (2004), Schmcja & Klessen (2004), Bonnell et al. (2006), and Bon- 
nell & Bate (2006). These and a number of other studies have led to two 
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suggestions in which a self-gravitating core could possibly evolve, viz. via 
gravitational collapse, or competitive accretion. While the former scenario 
envisages a self-gravitating core breaking up into a number of protostars, the 
latter invokes a baby protostar, typically ^0.5 Mq, that grows by accreting 
gas from a common pool while competing with its siblings to acquire mass. 
The scheme adopted in this work is akin to the latter scenario, where the 
mass of a protostellar baby, Minu, is taken to be ~ 0.05 M©, the Jeans mass 
at ~7 K, the lowest choice of prestellar core temperatures in our sample. 

The sound speed within a core, alf^ ~ al+a'^, where = hs^r^^ ff^ jg ^j-^g 
mean mass of gas molecules, for simplicity we assume the gas to be composed 
of molecular hydrogen and helium only; o" = for non-turbulent cores. We 
wish to emphasise that only those cores that are gravitationally bound are 
allowed to fragment, unbound cores are considered still-born and therefore 
disqualified for further calculations. We, however, do not accommodate the 
converse possibility of a super-Jeans core remaining starless although such 
cases have been reported (e.g. Sadavoy et. al. 2010). 

Protostellar fragments in either scenarios, gravitational fragmentation or 
competitive accretion, grow by accreting gas, calculations by Larson-Penston 
and Shu arrive at an accretion rate - 

^3 

Mace OC -g^, (4) 

which is quite different from the Bondi-Hoyle(BH) accretion rate. 



M. 



acc(BH) ~ T- 

^3 



frag 



(Bondi 1952), Vg being the terminal velocity of gas of average density, p, 
accreted by the fragment; Mfrag is the mass of the fragment. Schmeja & 



Klessen (2004), with their suite of numerical simulations suggested an em- 
pirical time-dependant rate of accretion 

log(M„ec) = log(Mo)-^t exp(-i//3), (5) 

for fitting parameters Mq and /3. While the BH- accretion rate provides an 
upper limit, that due to Larson- Penston, Hunter or Shu, is constant in time. 
In case of a self-gravitating core, the mass, M(r), within a radius r is 

/r'=rcore r , Jlj-J. 

where the integral in the second equality on the right hand side above ac- 
counts for any modulation due to protostellar feedback. 

The foregoing discussion seeks to underline the difficulties in estimat- 
ing the rate at which a protostellar fragment accretes gas. In the simplest 
case, that of BH accretion, Macc{BH) ~ 10"^ yr~^, provides the upper 
limit, while that due to Larson- Penston and Shu is ~ 10~^ yr~^, which 
probably, is the lower limit. Observationally determined values of Mace are 
typically an order of magnitude higher in comparison to the latter (e.g. Bon- 
temps et al. 1996). In cases 1, 2, 3 and 4a, we therefore adopt a constant rate 
of accretion. Mace ~ 10"^ Mq yr~^, corrected for the mass lost via feedback, 
as discussed in §2.1 below. 

The mass of a protostar, Mfrag, is then 

dM(t) , 

Mfrag ~ Mi^U + —^^t, (6) 

where dM(t) /dt — Mace — Mioss, Mioss being the rate of mass loss due to pro- 
tostellar feedback. Each fragment, depending on its mass, Mfrag, is assigned 
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a temperature, Tj^ag-, obtained from the well-known pre-main sequence evo- 
lution tracks. In the present exercise, Tj^ag varies between ~(1500-6500) K, 
for Mfrag between ~(0.08-2.00) Mq. Feedback from a protostar ejects some 
mass from the natal core, and heats up the rest of it which is then tested for 
possible formation of the next protostar. The process is continued till the 
core is exhausted whence star-formation is quenched. 

2.1. Feedback and cooling 

Young protostars drive energetic winds, and often, are associated with 
more energetic, and collimated molecular jets which inject considerable en- 
ergy in star-forming regions and supposedly regulate star-formation episodes 
(e.g. Lada 1985). While energetic jets are likely to puncture the cocoon en- 
closing a protostar to deposit the associated momentum outside the parent 
core, the weaker protostellar wind, on the other hand, may have a significant 
interaction with collapsing core and heat it (e.g. Wilkin & Stabler 1998). 
Below we estimate the rate at which protostellar wind drives out gas from 
the collapsing core with the help of a model suggested by Wilkin & Stabler 
(1998). The model envisages feedback in the form a shell driven by the wind. 
The luminosity of the protostellar fragment as it accretes gas from the natal 
core is. 



^this equation follows from a trivial re-arrangement of the expression for energy of the 
driving the wind, 2Eacc ^ V^Mfrag ^ J Maccdt. 








acc 



The velocity, V^, of the protostellar wind is 




(7) 



10 



R = R(t = t growth), being the radius of the shelL If M/oss is the rate of mass 
loss due to the wind, then a key parameter, the mass transport rate, deter- 
mines if the protostellar wind could possibly escape from a self-gravitating 
core; it is the ratio of Mioss against the rate of protostellar infall, Mj, 

„ = ^.0.5(^-)' (8) 

where V// is the free-fall velocity!^; Mi = dM{t)/dt, and in the absence of 
feedback, Mj is as defined by Eqn. (4) above with 0.975 being the constant 
of proportionality (Shu 1977). We use Eqns. (6) and (8) to calculate the 
mass of a protostar. 

Cooling The temperature of the envelope of gas cocooning the young pro- 
tostar, and heated by it is Tgas ~ Tfrage^^ , t being the optical depth within 
the core and r ~ Nh2S- For a core predominantly composed of molecular 
hydrogen, its column density, Nh2, may be estimated using an empirical 
relation, 

~ 4.41 X 10^1 ( It^) cm-2 (Tielens 2005), 
VK • km/s/ 

and a typical absorption cross-section, s ~ 10~^^ cm~^. The heated gas 
within a self-gravitating core is allowed to cool via stimulated dust emission 
at infrared wavelengths. Although dust is the dominant coolant at typi- 
cal protostellar densities, collisions between gas molecules and dust grains 
also contribute towards cooling, however, the contribution due to the latter, 

^We borrow Eqns. (8) and (11) from Stabler & Wilkin (1998); respectively, M, ~ 
AttRVw<J shell, and a shell ~ ^TTMTT' Eliminating (t shell from the two equations for a spher- 
ically symmetric wind, we arrive at Eqn. (8). 
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defined by the cooling function, 



A, 



■gas— dust OC 



/T \ 0.5 

/ -Lgas \ ,rp _rp \ 

^ y y-'- gas -'-dust) 



(9) 



(Goldsmith 2001), is less significant compared to the former. Cooling due to 
dust is quantified by the corresponding cooling function. 



(Goldsmith 2001). 

The new temperature of the core, T„e«,, that has radiated excess of its 
energy via emission from heated dust grains is. 



The residual gas in a core having acquired this new temperature, the revised 
Jeans mass, Mjeans{T — Tnew), was calculated for it. The calculations were 
terminated in case Mj^ansiT = T„ew) exceeded the available mass within the 
self-gravitating core; feedback from the relatively massive protostars may, 
however, disrupt the natal core. Heating of prestellar cores due to cosmic 
rays has been neglected in these calculations, for their effect is hkely to be 
restricted to the periphery of a typical core. 

3. Results 

The relationship between the CMF and the IMF has previously been dis- 
cussed on purely statistical grounds by for instance, Goodwin et al. (2008), 
where the authors commencing from the CMF for the OMC suggested by 
Nutter & Ward-Thompson (2007), derived the stellar IMF by randomly sam- 
pling the CMF with different choices of the core-to-star efficiency, SFE. The 
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Figure 2: Starting with an initial population of cores distributed according to the power- 
law, see Eqn. (1) above, the derived mass distribution for fragments has been plotted for 
three realisations; see text for description. The universal IMF, and a lognormal fit to the 
derived distribution has also been plotted for comparison purposes. 




Figure 3: The efficiency of star- formation plotted collectively for all cores corresponding 
to the ten choices of temperature, Tcore- The plot on the left is for non-turbulent cores 
while that on the right is for those supported by transonic turbulence. The addition of 
turbulence has little impact on the SFE, as is amply demonstrated by the plots. See text 
for description. 
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IMF so derived, for an SFE~ 0.15, i.e. a 15% core-to-star conversion ef- 
ficiency, matched the universal IMF suggested by Kroupa (2002). In the 
present work, however, having invoked physical details of a typical star- 
formation episode, modulo the assumptions listed above, we have adopted 
a different route in our attempt to deduce a distribution of stellar masses. 
However, our treatment of the problem is different from the hierarchical frag- 
mentation scheme suggested by Larson (1973). While the hierarchical frag- 
mentation scheme produced a distribution similar to the canonical IMF, the 
type of fragmentation envisaged in it lacks supporting numerical evidence. 
We shall now proceed to discuss individual realisation in this work. 
Case 1: (Non-turbulent cores) We obtained a collective sample of ~166,000 
stars for various choices of core temperature, Tco^e, between 7K-30K. The dis- 
tribution of fragment masses for this case has been plotted as a histogram in 
red in Fig. 2. It is characterised by a knee at ~0.5 Mq, before steepening for 
masses rightward of the knee, acquiring a Salpeter-like nature. For fragments 
with Mfrag/MQ < 0.5, the distribution becomes rather shallow and peters off 
steeply into the dwarf regime. The mass of the smallest fragment obtained 
in this case is Mfrag/MQ ~ 0.08. While a demonstrative calculation such as 
the one proposed here is likely to introduce inaccuracies in the results, the 
distribution bears remarkable similarity to the universal IMF defined as. 




(12) 
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(Kroupa 2002), plotted using continuous bold, black segments. Also plotted, 
for comparative purposes is the lognormal IMF suggested by Chabrier (2003), 

iV(M,,.,).exp[ -"°«^7r'°^"°'1 . (13) 

for Mq ~ 0.08 M0, and 2(j^ ~ 1, that also fits the red histogram, the 
distribution of stellar masses for this case, quite well. 

Knowing the number of fragments produced by a core, the core-to-star 
efficiency, SFE, is calculated as- 

SFE = ^ ^^/'-"g . (14) 

The distribution of the SFE as a function of their masses and initial temper- 
ature, Tcare, IS shown in the left-hand panel of Fig. 3. The SFE, in general, 
lies between 5% to 15%, while a few cores that are relatively large convert 
up to ~25% of their gas into stellar objects. Small, warm cores, on the other 
hand remain sterile as expected. 

Case 2: (Turbulent cores) Additional turbulent support within the precol- 
lapse cores, which in this case is transonic, not only produces fewer fragments 
(~ 159,000), but also shifts the distribution of fragment masses rightward, 
towards a higher mass, as is evident from the blue histogram shown in Fig. 
2. Greater turbulent support raises the thermal Jeans mass, manifested by a 
rightward shift in the knee of the resultant stellar mass distribution that now 
appears at ~0.8 M©; it falls off rapidly for Mf^ag/MQ >0.8. The formation 
of dwarfs though, according to the simple arguments presented here, may 
not be significantly affected due to the introduction of an additional thermal 
component in the form of a transonic turbulent velocity field. The arguments, 
nevertheless bring us to the interesting question of the fate of protostellar 
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disks in turbulent ambiance. Observational evidence in this regard is scant 
for lack of sensitivity of devices available at present (James di Francesco 
2010; Jenny Hatchell 2010 - private communication). One may have to wait 
for the commissioning of ALMA that promises to offer higher sensitivity and 
resolution. Despite the additional turbulent support, the core-to-star conver- 
sion efficiency, however, does not show much change in comparison to that 
for case 1, apart from a slight increase in the number of cores either with a 
low SFE, i.e. between 5% - 10%, or those that remain sterile; see Fig. 3. 
Case 3: (non-turbulent cores + no cooling) Cores, in this case, were al- 
lowed to fragment without invoking the cooling mechanism described in §2.1, 
although feedback from the prostellar fragments was included and ended up 
with even fewer fragments, ~ 156, 000. The effect of neglecting emission from 
warm dust in a self-gravitating core is similar to that of adding a turbulent 
velocity field to cores, reported in the previous case. The resulting distribu- 
tion of fragment masses for this case is shown by the yellow histogram plotted 
in Fig. 2. As in the previous case, while the high-mass end, Mfrag/M^ > 
0.8, is little affected, its low-mass end shows a considerable rightward shift. 
There are no fragments less massive than ~0.1 Mq. 

Case 4a: (non-turbulent cores -|- cooling; uniformly distributed 
core masses) Starting with a sample of 20,000 cores with masses in the 
range (0.5,10)Mo, but now distributed uniformly, the stellar mass distribu- 
tion for about 250,000 fragments derived using thermodynamic details as in 
case 1, has been plotted in the top panel of Fig. 4. Overlaid on top of it is 



16 



the power-law, 



dN{Mfrag) OC Mj^l/Mfrag 



a =-0.7; 0.08 < 
a = -1.5; 0.08 < (^) < 0.5 (15) 
a = 2.35; 0.5 < (^) 



which is not too different from the universal IMF defined by Eqn. (12) above, 
making it a rather interesting case. The knee of this distribution, for instance, 
is still located at M/^a^/M© ~ 0.5. It therefore leads us to the conjecture 
that the original nature i.e. the shape of the distribution of prcstcllar core 
masses may not be crucial to the deduction of the stellar IMF, and rather, 
the physical processes leading to the birth of stars are hkely to determine the 
nature of stellar mass distribution. The often reported similarity between the 
CMF, and the IMF, is therefore likely to be a matter of mere coincidence or 
a result of the incompleteness of our surveys of young star-forming regions. 

Case 4b (Same as 4a, but higher rate of accretion) To test the role 
of physical processes involved, we repeated case (4a) with a higher rate of 
accretion, M^cc ~ 10~^ Mq yr~^. The resulting distribution of stellar masses 
is shown in the bottom panel of Fig. 4; note that the distribution is signif- 
icantly overpopulated in the range (O.2,1)M0. Thus, protostellar fragments 
accreting matter with greater efficiency tends to flatten the resulting distri- 
bution of stellar masses. The contrasting results from these two reahsations 
re-emphasise the critical role of the physical environment in which protostars 
form for the rate with which gas is accreted from their natal cores depends 
on, both, the prevalent physical conditions within the core, and the rate at 
which mass is blown away by feedback. 
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Figure 4: Top panel /The distribution of stellar masses derived in the case commencing 
with a uniform distribution of core masses. The power-law given by Eqn. (15) has been 
overlaid for comparison purposes. Bottom panel : The stellar mass distribution derived for 
conditions similar to those for the results of the realisation plotted in the panel above. The 
accretion rate in this case is an order of magnitude higher in comparison to the former. 
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4. Discussion 

Determining the origin of the stellar IMF is a problem fundamental in na- 
ture, with arguments, both, supporting and confuting the hypothesis tested 
in the present work. Literary evidence favouring the hypothesis is generally 
encouraged by the apparent similarity between the observationally deduced 
CMF for a number of star-forming regions, and the stellar IMF (e.g. Motte 
et al. 1998, Luhman et al. 2006, Ward- Thompson & Nutter 2007). On the 
other hand, Bonnell et al. (2006), for instance, on the basis of their nu- 
merical simulations of self-gravitating molecular clouds argued that physical 
processes leading to the interplay between self-gravity and the thermal pres- 
sure, and possibly the magnetic field (e.g. Kirk et al. 2009), could perchance, 
determine the nature of stellar mass distribution. In this work we examined 
the hypothesis in light of a simple semi-analytic scheme outlined in §2 above. 

The four cases discussed above seek to emphasise the importance of ther- 
modynamics and the rate of accretion in a self-gravitating core. While an 
accretion rate of ~ 10~^ Mq yr~^ in cases 1, 2, 3 and 4a led to a distribution 
that resembled the universal IMF, increasing it by an order of magnitude 
produced a flat distribution of stellar masses; see plot for case 4b in the 
bottom panel of Fig. 4. We therefore, as evidenced by results in case 4a, 
argue that a power-law CMF is not a necessary preclude to a stellar mass 
distribution that resembles the universal IMF. 

We note that the models in cases 1,2, and 4a which include contribution 
due to dust coohng and a relatively lower rate of accretion. Mace ~ 10~^ M© 
yr""^, irrespective of additional turbulent support, produce dwarfs as small 
as ~0.08 Mq, however, elimination of dust-cooling from calculations in case 
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3 shifts the stellar mass distribution towards a higher mass, with a knee 
at 0.8 M0 without any fragments smaller than ~0.1 Mq. We allude this 
deficiency to the upward revision of the initial protostcllar mass, Minu, due 
to warmer gas. Remarkably enough, starting from a uniform distribution of 
core masses in case 4a, but with thermal details identical to those in case 1, 
the stellar mass distribution derived here resembles the universal IMF which 
emphasises the relative importance of physical processes such as the rate of 
protostellar accretion over the initial distribution of core masses. Increasing 
Mace, as in case 4b, for instance, led to a considerable overpopulation in 
the range (O.2,1)M0, rendering the distribution flat. Having demonstrated 
that a power-law CMF of the type defined by Eqn. (1) above need not 
necessarily preclude a stellar mass distribution that resembles the universal 
IMF, the argument favouring a direct one-to-one relationship between the 
two is called into question; we do not envisage one on this line either. 

Cores, in general, appear to convert between 5% to 15% of their gas in 
to stellar fragments, while a few cores of intermediate masses may have a 
gas-to-star conversion efficiency of up to 25 % as can be seen from plots of 
the SFE in Figs. 3. The SFE for a sample of cores, as defined by Eqn. (14) 
above appears agnostic to changes in ambient conditions within putative star- 
forming cores, albeit the evolution of an individual core is governed critically 
by the prevalent ambient conditions. The SFE derived here is consistent 
with that reported for star-forming regions in the Serpens, Perseus, and the 
Ophiuchus molecular clouds (e.g. Enoch et al. 2008, and references therein). 
Despite the success, the scheme is too naive to encompass the complexity of 
the star-forming process. It may also suffer owing to the lack of a complete 
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understanding of the formation of dwarfs; possible explanations for which 
range from ejection at the embryonic stages (e.g Reipurth & Clarke 2001), 
or via fragmentation of protoplanetary disks (e.g. Goodwin & Whitworth 
2007, Stamatellos & Whitworth 2009). 

Finally, we must consider the possible effect of the timescale on which a 
core evolves, on the distribution of stellar masses. The problem was broached 
by Clark et al. (2007), the prestellar cores in the present calculations, how- 
ever, are immune to this problem as they all had roughly similar initial 
densities, and therefore evolve on a mutually comparable timescale. In view 
of the empirical nature of Larson's scaling relations and their application 
over a large spatial range, we wish to suggest that the time-scale problem as 
enunciated by the above authors may not really be a nemesis to the stellar 
IMF. 

5. Conclusions 

For two samples of prestellar cores, distributed according to : (A) a 
power-law distribution defined by Eqn. (1) above, and (B) uniformly, with 
masses in the range 0.5 M0 to 10 Mq, we have derived a distribution of 
stellar masses. We draw the following inferences : 

1. By invoking a simple cooling mechanism and protostellar feedback, 
a stellar mass distribution similar to the IMF, is recovered for either 
choices of the CMF, (A) or (B), with M^cc ~ 10~^ Mq yr^^ consistent 
with that inferred observationally; a higher rate of accretion, however, 
produces too many sub- Solar protostars as can be seen in the bottom 
panel of Fig. 4, and could thus, be a crucial parameter in determining 
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the final shape of stellar mass distribution. 

2. The distribution of stellar masses appears sensitive to the thermal prop- 
erties of self-gravitating cores, and the rate at which protostars ac- 
crete matter from the natal cores (Cases 1, 2 & 3 in §3), interestingly 
enough, it appears intransitive to the preceding distribution of core 
masses (Cases 4a & 4b). We therefore suggest, the apparent similarity 
between the observationally derived CMF, such as the one given by 
Eqn. (1), and the universal IMF is possibly a coincidence. 

3. The SFE, however, appears unaffected by the injection of a little tur- 
bulence in parent cores (case 2 in §3) , although the stellar mass distri- 
bution is shifted a little towards the higher mass end, and the shift is 
considerable when self-gravitating cores arc not allowed to cool. There 
is a deficiency of fragments less massive than ~0.1 M0 in the latter 
case. 

4. The prestellar cores in our test sample had roughly similar initial den- 
sity modulo variation by a factor of 2-3, and so, they evolved on a 
mutually comparable timcscalc. In view of the empirical nature of the 
scaling relations defined by Eqns. (2) and (3) above, it appears that 
despite spanning a wide range of masses, the prestellar cores may have 
approximately comparable initial densities, and therefore, evolve over 
similar timescales. The distribution of stellar masses derived here is 
therefore fairly robust, and inoculated against variations in the evo- 
lutionary timescales of the cores in the ensemble. In order to place 
this argument on a stronger footing, we must further investigate the 
formation of prestellar cores, and their congenital properties. 
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